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We suggest a method of fitting of a zero-range model of a tectonic plate under 
a boundary stress on the basis of comparison of the theoretical formulae 
for the corresponding eigenfunctions/eigenvalues with the results extraction 
^ under monitoring, in the remote zone, of non-random (regular) oscillations 

^ of the Earth with periods 0.2-6 hours, on the background seismic process, in 

case of low seismic activity. Observations of changes of the characteristics 
^ of the oscillations (frequency, amplitude and polarization) in course of time, 

G together with the theoretical analysis of the fitted model, would enable us to 

localize the stressed zone on the boundary of the plate and estimate the risk 
^ of a powerful earthquake at the zone. 
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O 1 Dynamics of the system of tectonic plates 

O and the motivation of the zero-range model. 
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The lithosphere of Earth consists of 14 tectonic plates which jigsaw fit each 
other. The plates are isolated from underlying solid structures within the 
Earth mantle by the low- viscosity layer of the asthenosphere which is formed, 
due to various conditions, including high pressure and temperature, in the 
interval of depth 100 -200 km. The plates move on the surface of Earth, due 
to convective flows and variations of the angular speed of Earth, gliding on 
the low-viscosity layer of asthenosphere and interacting with each other at 
some active boundary zones, see below. Fig. 1 below shows a complex form 



^ Fig.l is borrowed from the book of S. Aplonov, [3 , and is included in our text with 
permission of the Pubhshing House of the St. Petersburg University. 
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Figure 1: 

Boundaries of tectonic plates: 1- divergent boundaries (a - oceanic ridges, b- 
continental rifts), 2- transforming boundaries, 3-convergent boundaries (a- 
insular, b- active continental outskirts, c - collisions of plates). Directions 
and velocities of movement of plates ( cm/ year). 

of boundaries between tectonic plates and their movements in different direc- 
tions. The plates are enumerated in the following order : 1. South- American 
plate, 2. African plate , 3. Somali plate, 4. Indian and Australian plates, 
5. Pacific plate, 6. Nazca plate, 7. North-American plate, 8. Philippines 
plate, 9. Euro-Asian plate, 10. Antarctic plate, 11. Caribbean plate, 12. 
Cocos plate, 13. Arabian plate. 

In addition to the 14 large plates described above there are 38 smaller 
plates (Okhotsk, Amur, Yangtze, Okinawa, Sunda, Burma, Molucca Sea, 
Banda Sea, Timor, Birds Head, Maoke, Caroline, Mariana, North Bismarck, 
Manus, South Bismarck, Solomon Sea, Woodlark, New Hebrides, Conway 
Reef, Balmoral Reef, Futuna, Niuafo'ou, Tonga, Kermadec, Rivera, Gala- 
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pagos, Easter, Juan Fernandez, Panama, North Andes, Altiplano, Shetland, 
Scotia, Sandwich, Aegean Sea, Anatoha, Somaha), for a total of 52 plates. 

It is commonly accepted that the tectonic plates are relatively thin elas- 
tic structures, approximately 100 km. thick, with linear size from 1000 km 
to several thousand km. The material of the plates, at the depth 100 km., 
has typical Young's modulus 17.28 x 10^° kg mT^ sec'"^ , density 3380 kgm^^ 
and Poisson coefficient 0.28. The velocity of the longitudinal waves in these 
materials is approximately 8000 m sec~^, and the velocity of the transversal 
(fiexural) eigen-waves depends on the eigenvalues and varies, depending on 
the type of the wave, on a wide range around 4500 m sec~^, see the formula 
^ below. Because of non accurate matching of the boundaries of the neigh- 
boring plates, the zones of direct contact of the plates are typically small- 
about 100 km.- compared with the linear size of the plates. Remaining inter- 
plate space is filled with loose materials, which are not able to accumulate 
any essential amount of elastic energy caused by the deformation. These 
materials can damp oscillations with short periods as 20 min. - 1 hour. 
Damping of acoustic waves by loose materials was discussed in [I9l [HI |35] . 
Since tectonic plates are relatively thin, a major part of their kinetic energy 
is stored in the form of oscillatory fiexural modes. The underlying layer of 
the asthenosphere makes fiexural oscillation of plate possible, but also helps 
damping of the fiexural waves with short periods, due to non-zero viscosity. 
Based on above data we conjecture that the fiexural modes in different tec- 
tonic plates, in a certain range of periods, supposedly between 0.2h — Ih, see 
below the estimation of periods of eigen-modes of a model rectangular plate, 
are elastically disconnected from each other. Thus we expect that these fiex- 
ural modes characterize elastic properties of the plates, but not the global 
elastic properties of the Earth's crust. 

The convective fiows in the asthenosphere and variations of the angular 
speed of Earth, thanks to long-time variations of the moment of inertia of 
Earth, may cause collisions of neighboring plates. In presence of the liquid 
friction on the underlying layer of the asthenosphere, variations of the angular 
speed of earth cause mutual displacements of the neighboring plates, because 
small plates react immediately on the variations of the angular speed, and 
larger plates lag behind. These displacements cause collisions of the plates 
in the active zones, where the plates directly contact each other. Generally, 
when the angular speed of Earth decreases, with growing of the moment of 
inertia due to displacement of the center of gravity of Earth, the collisions 
may occur on the eastern boundaries of the major plates contacting smaller 
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plates, observe, for instance, the contact of the Euro-Asia plate and the 
Philippines plate on Fig. 1. The stress caused by the collision may be either 
discharged due to forming cracks in the plates, splitting the active zone into 
independently moving fragments, or, being applied for an extended period 
of time, may cause accumulation of a considerable amount of (potential) 
elastic energy in the active zones of contacts, in form of elastic deformation 
of the stressed plates. This energy may be eventually discharged in form of 
a powerful earthquake. 

Accumulation of the elastic energy, due to standard variational princi- 
ple [7j, causes the increment of eigenfrequencies of the tectonic plates. In 
[271 1321 E3] oscillatory processes, with similar spectral characteristics, were 
observed in mutually remote zones of the Euro-Asia plate. The authors of 
[27] suggested calling the processes seismo-gravitational oscillations of the 
Earth ( SGO). Recent analysis, based on data of the international GEO- 
SCOPE network, revealed the existence of global oscillations with periods 
3.97 h, 3.42 h, 1.03 h and 0.98 h. SGO with smaller periods do not have 
global character. They are usually observed in certain plates. For measure- 
ments of SGO, special devices are used ("vertical pendulum"), with high 
sensitivity to the variations of amplitudes and frequencies of SGO within the 
interval of periods 0.2 - 2 hours. Modification of the registering channel of 
the device allowed us to extend the interval to 0.2-6 hours. Results of analy- 
sis of measurements of SGO with relatively short periods from this range are 
presented in [30]. The data on SGO with longer periods and characteristics 
of modern 3-channel seismographs can be found in [29j. Observations with 
these devices reveal a wide spectrum of SGO. Some of the measured fre- 
quencies of SGO coincide with short-time variations of the angular speed of 
Earth, which were estimated based on astronomical observations, see the [2B] ■ 
All experimental data confirm the presence of active energy in the system 
of tectonic plates, relevant to the inhomogeneity of the lithosphere. Direct 
calculations of eigenfrequencies of the fiexural modes of a rectangular thin 
plate f2 = c?! X (^2, (ii = 4 000 km, = 8 000 km, 200 km thick, were done, 
see [21] , based on bi-harmonic model, under natural assumptions concerning 
the density of the material is 3380 kg/m^, the Young's modulus is 17.28 x 
10^° kgm~^ sec'"^, and the Poisson coefficient is equal to 0.28. Though the 
Young's modulus of a tectonic plate varies on a wide range of values, it is 
possible to provide a rough estimation, see below, of the periods of SGO and 
the length of the corresponding running waves of a tectonic plate, by comput- 
ing these parameters for a model rectangular plate with elastic parameters 
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Figure 2: 

Computed eigenfrequencies of the rectangular plate 4 000 km x 8 000 km 
for the modes f(l, N) 

equal to those of the tectonic plate on the half depth, 100 km. Indeed, the 
corresponding dynamical equation for the transversal (vertical) displacement 
/is 



du 



an 



dAu 
dn 



= 

an 



In PI] the simplest Neumann boundary conditions 

are imposed on the boundary. This allows us to solve the dynamical equa- 
tion by Fourier method, via separation of variables. The eigenfunctions of 
the bi-harmonic operator coincide with the eigenfunction of the Neumann 
Laplacian, /^^ /j = cos ^^^^p- cos but the eigenvalues of the Neumann 

bi-harmonic operator 

are squares of the corresponding eigenvalues of the Neumann Laplacian. Tak- 
ing into account that 

0.167 10^0= [0.64 xlO^]'=:a\ 



12p(l - a 
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we find the periods of flexural eigen-ocsillations from the formula: 



or 

271- 

^ii,h = Tf; — = >^th =■ ^Li2' (2) 

ll,l2 

where k plays the role of the corresponding momentum. The speed of the 
transversal ( flexural) waves is calculated as 

= = 2 a ki,,i,. (3) 

Then for the above model plate the periods of flexural oscillations are defined 
by the formula Thours — 2.63 [Alf + /|]~^ hours, and the velocity of the cor- 
responding flexural waves are calculated as vj ^, m sec . 
In particular, the period of the flexural oscillation /i^2 and the speed of the 
corresponding flexural wave v{^2 calculated as 

ri,2 = 0.33 hours, v^^ = 9520 m sec~\ (4) 

and less for longer periods. The corresponding space-temporal oscillatory 
mode is 

TTXi 7r2x2 

COS uii COS cos , 

4000 8000 ' 

with Xi, X2 measured in kilometers. This oscillatory mode can be represented 
as a linear combination of running flexural waves 



cos 



n Xi 71X2 
Uiot ± ± 

4000 4000 



(5) 



The corresponding wavelength is estimated by the minimum of v^A? + Ai, 
where Ai, A2 correspond to the shift of the running wave by the corresponding 
period in time. For instance 

TT Ai + vrAs 
4000 

gives an estimation of the minimal wavelength as -^/A^ + A^ > 5656 km, 
which is much more than the diameter of the active zone, 5656 km >> 100 



6 



km. Hence the zero-range model can be used, under above assumption, for 
the model stressed plate. 

It appeared that periods of the eigen-modes — /s^s sit in the inter- 
val 0.21 h - 5 h and their total number and distribution looks similar, see 
Fig. |2| to SGO described in the paper [271 1321 E3], despite the trivial plane 
rectangular geometry of the plate and trivial Neumann boundary conditions. 

Examination of the results observed in [271 ESI EH] and theoretically ob- 
tained in [21] resulted in the conjecture that transversal SGO can be inter- 
preted as flexural eigen-modes of the relatively thin tectonic plates. We hope 
that the zero-range model may be also used in this periods range for real 
tectonic plates as well. 

In the next section we model dynamics of tectonic plates based on the 
bi-harmonic boundary problem with "natural" boundary conditions. 
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Figure 3: Dynamics of seismo-gravitational oscillations observed in St. Pe- 
tersburg, 22-29 March 2000, see the text below. 

Variation of the character of SGO is seen from comparison of the graphs 
of the amplitudes of SGO on 48-hours intervals of time separated by the 
vertical dotted lines. The graphs on these intervals show the reaction of the 
device on the oscillations of the ground. 

A typical phenomenon of intense seismo-gravitational pulsation (SGP) 
was registered on the vertical component Z on the initial interval marked by 
0. This phenomenon was also noticed in the earlier paper [13]. See more 
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about SPG in the further text, after the figure |4j 

The components Z,EW and NS are obtained after filtration. The graph 
on the interval 1 shows the reaction of the filter on the maximal phase of the 
earthquake 28 March 2000 in Japan. The non-filtered graph Z represents the 
reaction of the base of the vertical seismograph on the transversal oscillations 
of the ground. 

Spectral-time cards, see below fig. |4] represent the effect of growing of 
frequencies of SGO in visual form. The data of the observations were filtered 
by the band filter with the band strip [60 min.-300 min.]. We used the 
gliding time window length 700 min., and step- wise shifts with the 5 minutes 
steps. The interval of frequencies, measured in micro-Hertz, was chosen as 
[70mcHz,250'mcHz] =: [Fmin, Fmax]- For resolution O.OSAmcHz the steps 
were halved to guarantee better smoothness of the spectral function. 

Effect of growing of the frequency of SGO was detected only on spectral- 
time cards (ST-cards) of the vertical component of SGO, see below fig. |4] 
The frequencies are marked on the horizontal axes, time in hours - on the ver- 
tical axes, for the intervals 0,1,2 respectively. The level of spectral amplitudes 
is represented by the variation of the color. Oscillations with large spectral 
amplitudes, higher than the average level on the spectrum, are marked by 
grey color. The maximal amplitudes are marked by white color. The values 
of these amplitudes on the interval are 3.5 times greater that the average 
amplitude, ( frequency about 200 mcHz). On the interval 1 they exceed the 
average amplitude in 1.8 times (the frequency about 200 mcHz). On the in- 
terval 2 they exceed the average amplitude in 3 times ( the frequency about 
160 mcHz). For SGO with frequencies 90-110 mcHz and 170- 190 mcHz the 
amplitudes are 1.75 and 2.75 is times larger than the average amplitude. All 
three parts 0,1,2 of the fig. |4] reveal two patterns of inclined gray domains 
( left-down to right-up and left-up to right-down) which correspond to SGO 
with growing and decreasing frequencies, correspondingly. This patterns pro- 
vide an evidence of increment of the stored elastic energy in the system and 
discharging the stored elastic energy, respectively in form of oscillation modes 
with certain frequencies. 

Growing of the frequency on the interval 2: 

1. The frequency is growing from fmin = SOmcHz to 123mcHz during 
the period of 15 hours. 

2. The frequency is growing from fmin = 130mcHz to 185mcHz during 
the period of 28.5 hours. 

3. The frequency is growing from fmin = l70mcHz to 200mcHz during 
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the period of 8 hours. 

Growing of the frequency on the interval 

1. The frequency is growing from fmin = 165mcHz to 205171011 z during 
the period of 33 hours. 

2. The frequency is growing from = 12hmcHz to IQthmcHz during 
the period of 15 hours. 

Growing and decreasing of frequencies of SGO are easily noticeable also 
on the interval 1. 

Growing of the frequency is characterized by the ratio where At/ is 
the increment of the frequency and r is the corresponding time interval, when 
the growing was observed. 




80 100 120 140 160 180 200 Wi 



Figure 4: Time-frequency analysis of the vertical component of SGO 

Other domains, where the frequency of the modes decrease with time 
growing, may be interpreted as an evidence of local relaxation of the stress, 
probably caused by local destruction of the plate (forming cracks). 

We conjecture that the extended growth of frequencies of the SGO may 
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be considered as a precursor of strong earthquakes. Our ability to extract 
useful information from the observations of frequencies and the shape of 
SGO is limited by our understanding of the mechanism of the variation of 
the frequency and the shape of SGO modes arising from the boundary stress 
on the tectonic plates. 

There were also other observations of the shift of frequencies of selected 
modes of SGO. Note, that in numerous observations on SGO in Leningrad 
(St. Petersburg) intense short pulses were also recorded. They are consti- 
tuted by several sinusoidal harmonics with periods from 30 minutes to 1 
hour, and the total duration of the process 6- 10 hours. In several cases they 
were also followed by powerful earthquakes - in 2-4 days. This process was 
noticed first in [13] and was given the name of "seismo-gravitational pulsa- 
tions", SGP. Statistical analysis was done based on the data of the 6 months 
monitoring of SGO in St-Petersburg. This analysis confirmed that the con- 
nection between SGP and the subsequent strong earthquake is not random, 
with probability 95 %. 

Essential information on SGO is obtained from the observation of vari- 
ations of the frequencies and the shape of the corresponding modes in the 
remote zone. The typical size of the zone of contact is negligible - "point- 
wise"- when compared with the wavelength, which generally can be estimated 
as [A^]"^^^ see (Q). 

Although the boundaries of real tectonic plate are not smooth, however 
details of the local geometry of the boundaries may be neglected compared 
with the wavelength of typical waves on the plates. Then for SGO with 
periods 20 min. - 1 hour we may base the theoretical analysis of SGO on 
the bi-harmonic model for the relatively thin plate with natural boundary 
conditions. In this paper we also neglect presence of the liquid layer of the as- 
thenosphere underlying the plate, thus reducing the problem to construction 
of self-adjoint extensions of the 2-d bi-harmonic operator. 

In [31] we use the fact that characteristics of the seismo-gravitational 
oscillations in remote zone depend on a small number of basic parameters. 
On one hand, his fact is can be interpreted in a spirit of the Saint- Venant 
principle |^ On the other hand we were able to interpret this fact in the 
spirit of operator extensions. Based on this observation we developed in [31] 
a preliminary version of the above arguments and suggested to use a solvable 



^We are grateful to Doctor Colin Fox for inspiring discussion concerning the Saint- 
Venant principle 
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zero-range model of the tectonic plates under a point-wise boundary stress, 
caused by the colhsion of plates. The role of the parameter of the model 
was played by some 3x3 real Hcrmitian matrix M, see also next section. 
Thus the number of the Saint- Venant parameters for the boundary stressed 
relatively thin tectonic plate is 6. We presume, that the matrix M defines 
the type of the stress, depending on mutual positions of the contacting plates 
at the active zone. There is a good reason to call M the Saint- Venant ma- 
trix and the number of the Saint-Venant parameters (6) - the Saint-Venant 
number. We assume that mutual positions of the plates in active zone, and 
hence the matrix M remains essentially unchanged during extended period. 
Thus the matrix M characterizes the type of contact of tectonic plates in 
this location and remains the same for all earthquakes arising from the given 
active zone. In this paper, based on the compensation of singularities in the 
fundamental Krein formula, we suggest an explicit formula for the perturbed 
eigenfunctions of the plate. Comparison of the calculated eigenfunction with 
the data of the instrumental observations in remote zone, permits, in prin- 
ciple, to fit the model, that is, to find the basic matrix M. Once fitted, the 
constructed model would allow to calculate in explicit form the increment 
of the eigenvalues and the variations of the shape of the eigenfunctions of 
the stressed (perturbed) plate depending on the type and the magnitude 
of the local stress. We assume that the shifts of the eigenfrequencies and 
the changes of the shape of the eigen-modes can be measured for each active 
zone. With numerous active stations in the GEOSCOPE and IRIS networks, 
fitting of the proposed model can be done, eventually, for all active zones on 
the boundary of each major tectonic plate. 

Generically earthquake hits only one active zone at a time. Then com- 
paring the observed shift of the eigenfrequencies and the variations of the 
shape of the mode in the remote zone, where the Saint-Venant principle is 
applicable, with the results of computing based on the model, we will be 
able to locahze the excited active zone based on the observed changes of the 
eigen-frequencies and the shapes (amplitude, polarization) of the flexural 
eigen-modes. 

Mathematically the zero-range model of the isolated point-wise stressed 
tectonic plate and a similar zero-range model of the point-wise stressed tec- 
tonic plate submerged into environment formed by other plates, intermediate 
layers and asthenosphere, differ by the type of basic equations, but have a 
lot in common. In particular, the number of free parameters (6) for the 
point-wise active zone, in bi-harmonic model and Lame model is the same. 
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We interpreted these parameters in the spirit of the Saint- Venant principle 
as essential parameters describing the shape of the wave-process in the re- 
mote zone, see [31]. In this paper we obtain the first order approximation for 
the perturbed eigenfrequencies and eigenfunctions based on a modified Krein 
formula, see [IHl H], for the point-wise stressed thin plate described by the 
bi-harmonic equation. The explicit representation of the perturbed eigen- 
modes permits to fit the zero-range model suggested in [31] based on results 
of instrumental measurements of SGO. We also conjecture that the properly 
fitted solvable model of the stressed tectonic plate may help to enlighten the 
nature of the pulsations. 

2 Zero-range model of the point- wise bound- 
ary stress 

We may base our approach on the standard mathematical model of the tec- 
tonic plate in form of a thin elastic plate, thickness h, with free edge, or on 
the 3-d Lame equations for displacements. We consider both options, de- 
scribing in the next two subsections specific details of both models. Then we 
develop the common part of the theory, for both models simultaneously. 

2.1 Thin plate model for the isolated tectonic 
plate under the localized boundary stress 

Denoting by D the "the bending stiffness " , connected with Young modulus 
E, the thickness h of the plate and the Poisson coefficient a by the for- 
mula D = Eh'^[12{l — cr^)]^^, we represent, following [121, corresponding 
dynamical equation for the normal displacement u as 

We consider time-periodic solutions u{ujs, x)of the equation and separate the 
time, thereby reducing the dynamical problem to the spectral problem with 
the spectral parameter 

A = pD-^ (6) 
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for a bi-harmonic operator on a compact 2-d domain f2 - the tectonic plate 
- with a smooth boundary dQ: 



Au = A\ = Xu,ue W^{Q), 



and free boundary condition involving the tangential and normal derivatives 
of the displacement u and the tension Am: 



dAu 
dn 



+ (1-^) 



^3 



U 



dndf^ 



an 



0, 



Am - (1 - cr 



an 



0. 



(7) 



Here n, t are the normal and the tangent directions on the boundary. 

The bi-harmonic operator A is selfadjoint in the Hilbert space L2{^1) := 
7i. The eigenfunctions of A are smooth and they form an orthogonal basis 
in L2{fl) = Ti. We consider the restriction Aq of A onto D{Aq) constituted 
by all smooth functions vanishing near the boundary point a G dVL. The 
restriction is symmetric, but it is not selfadjoint, because the range of it 
[A — XI) D{Aq), for complex A has a nontrivial complement Nx which is a 
linear hull of the Green function G{x,a,X) := go{x,X) and its tangential 
derivatives ^'^^^'^'^^ ■= g^(^x,X), ^ '^g^f'''^^ := g2{x,X) of the first and second 
order, at the point a. The orthogonal complement A^a of the range is called 
the "deficiency subspace" , and elements of it - "deficiency elements" : 



Nx:=ne {Ao - XI)Do = \/ ^,(*, A). 



s=0 



The deficiency elements have, at the boundary point a, singularities of dif- 
ferent types ( see for instance [H], where much more general problem is 
considered) : 

go{n, t) ^ (n^ + 1^) \^(^n^ + 1'^)^ g^(^n, t) ^ t \n{n^ + 1^), g^{n, t) ^ \n{n'^ + f), 

hence they are linearly independent and form a basis in the deficiency sub- 
space. The deficiency subspace at the spectral point A is 



ne{Ao-XI)Do = \/ gs*,X. 



s=0 
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The dimensions of the deficiency subspaces (3,3) constitute the "deficiency 
index". Hereafter we select X = i and attempt to construct a self-adjoint 
extension of Ao, which will play a role of a zero-range model of the tectonic 
plate under the boundary strain. 

Note that Lame model the deficiency index is also (3,3), on a smooth 
boundary. The role of deficiency elements is played by the columns of the 
Green matrix. The boundary of the tectonic plates may be assumed smooth 
for the long waves (small A), since the integral shape of the solutions of the 
differential equations with small A is not affected by the details of the local 
geometry. 



2.2 Construction of the self-adjoint extension 

Extend Aq from Dq onto D{Aq) = Dq + Ni + N^i as an "adjoint operator" 
Aq by setting (Aq ± il) g = for g E N^i. This operator not selfadjoint, 
and it is not even symmetric, so that the boundary form 



{A+u, v) - {u, A+v) =Jiu,v) 



(8) 



does not vanish, generally, for u,v E D{Aq). One can rewrite (|8j) in more 
convenient form with using new symplectic coordinates with respect of a new 
basis in 



w: 



_ 1 

" 2 
1 

2i 



9s 



A + il 
A-il 



9s 



A 



9s 



A + il 
A-il 



9s 



A-il 
I 

~A-iI' 



Since A^gs 



i9s 



0, K- 2/]^^, = Owe have. 



(9) 



Following [21] we will use the representation of elements from the domain of 
the adjoint operator, by the expansion on the new basis: 



s s s 

A 



9s 



Uq + 



A-il 



A-il 



(10) 
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Note that due to Q 



A — il A — il A — il A — il 

Note that the boundary form 

(A+M, t;) - (m, A» ■.= J{u,v) 

of elements u,v, 

u = u^ + j^jil - -J^C- ■■= uo + n^ uo e DiAo) G iV, (11) 



A~ir+ A- if 
depends only on components n", ri" of -u, f in the defect A^. Then the bound- 
ary form is represented as: 

{A+u,v) - {u.Atv) := J{u,v) = {Ct,C) - {C,a) (12) 

with Euclidean dot-product for vectors C,± G Ni. Note that the representa- 
tion of the boundary form in terms of abstract boundary values C,± contains 
only integral characteristics of the elements from the domain of the operators 
considered, and hence it is stable with respect of minor local perturbations 
of geometry of the plates. This enables us to substitute, for practical calcu- 
lations, the real irregular boundaries of the plates by the smoothed bound- 
aries, obtained via elimination of minor geometrical details, compared with 
the length of standing waves, circa 5000 km, of SGO with periods in the 
essential gauge 20 min - 1 hour. 

The boundary form vanishes on the Lagrangian plane defined in D{Aq) 
defined by the "boundary condition" with an Hermitian operator M : ^ 

: 

ef = Me_. (13) 

This boundary condition defines a self-adjoint operator Am as a restriction 
of onto the Lagrangian plane Tm G D{Aq) defined by the boundary 
condition (13). The resolvent of Am defined by the boundary conditions is 
represented, at regular points of Am, by the Krein formula, see [21 l2l] : 

where P is an orthogonal projection onto N^. 
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2.3 Compensation of singularities in Krein formula 
and calculation of the perturbed spectral data 

Singularities of the resolvent {Am — A/)^^ coincide with the spectrum of 



Am- But both terms in the right side of (14) also have singularities on the 
spectrum of the non-perturbed operator A. The singularities of the first 
and second term the eigenvalues of A compensate each other. We are able 
to derive this statement via straightforward calculation, in classical Krein- 
Birman-Schwinger formula and in the corresponding construction for quan- 
tum networks, see p^fT]. In the course of the calculation of the compensation 
of singularities we can recover both the eigenvalues of the perturbed opera- 
tor Am and the corresponding eigenfunctions, see [18]. Note that a similar 
statement, as a lemma on compensation of singularities of the corresponding 
Weyl-Titchmarsh function, was discovered in [B] for 1-d solvable model of 
the quantum network in form of a quantum graph. Later, in [16] and in 
|15j . similar statements were proven for Dirichlet-to- Neumann maps of quan- 
tum networks. We formulate here this statement for the resolvent of the 
selfadjoint extension based on ideas proposed in [T7] . 

We will observe the effect of compensation of singularities on a certain 
spectral interval Aq = [Aq — 5, Aq + 5], centered at the resonance eigenvalue 
Ao of the non-perturbed plate, assuming that the perturbation defined by the 
matrix M is relatively small, in a certain sense, se below. 

Assuming that there is a single eigenvalue Aq of A on the interval Aq, 
with the eigenfunction ipo, we use the following representations, separating 
the polar terms from smooth operator functions Ki, K^i, K on Aq 

^ + M ^ -N^oK^ ^ ^ 

a3a7 = ^^° + ^)a^ + ^- 
a3a7 = ^^°-^)a^ + ^- 

P^PHl + Ag)^^^ + i^(A), (15) 

with a smooth matrix-function K = Kq + o(|A — Ao|), with Kq = K{\q) and 
II o(|A-Ao|) ||<Co<5. 

Definition 2.1 We say that the matrix M is relatively small, if [I+K{X)M]^^ 
exists and is bounded on Aq. 
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This condition is obviously fulfilled if 



\I + K{Xq)M]-' II Cq5 « 1 



(16) 



To calculate the second term in the right side of the Krein formula (14) we 



have to compute the inverse of the denominator, that is to solve the equation 

" = 9 (17) 



Though the standard analytic perturbation technique is still not applicable to 



this equation under the above conditions 2.1 or rtl6h, we are able to construct 



the inverse based on finite dimensionality (one-dimensionality) of the polar 
term. 



u 



[I + KM]-^g-{l + \ 



[I + KM]-^Pipo){PipoM[I + KM]-'g) 
\o-\+{P^oM[I + KM]'^ipo) ' 



Based on (18) we are able, see [HI [T] to observe the compensation of singu- 



larities in the above Krein formula (14) and calculate the polar term of the 



resolvent at the single eigenvalue of the operator Am on the interval Aq: 



Theorem 2.1 // the perturbation is relatively small, as required in (2.1), 



then there exist a single eigenvalue Xm of the perturbed operator Am on the 
interval Aq which is found as a zero of the denominator in [Tl 



(19) 



Ao - A + (1 + Xl){PipoM[I + KM]-'ifo) := dM(A), dM(AM) 
and the corresponding eigenfunction 

= ^o- (Ao - i)K_,M[I + KMY^Pip^, 



(20) 



computed at the zero Xm- The polar term of the resolvent of the perturbed 
operator at the eigenvalue is represented as: 

ipo - (Aq - t)K_iM[I + KM]-'Pipo){ipo - (Aq - t)K_,M[I + KM]-'Pipo 

dM{X) 

Proof of this statement can be obtained similarly to the corresponding 
statement in [IB], where the case of several eigenvalues of the non-perturbed 
operator on an essential spectral interval was discussed. We consider here the 
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simplest case, when only one eigenvalue of the unperturbed plate is present 
on an essential spectral interval. If the perturbation is small as required by 



the condition (16), then the approximate eigenvalue and the corresponding 
approximate eigenfunction of Am can be obtained via replacement K, K±i in 
( fl^ byi^'(Ao),/r±.(Ao): 



Am ^ Ao + (1 + \l){PipoM[I + K{\,)M]-^K{\,)^o), 

^^0- (Ao - t)K^ii\o)M[I + K{\o)M]~'Pipo. (21) 

Remark Analysis of the multi-point boundary condition which corresponds 
to several stresses applied at the points ai, 02, . . . ctm on the boundary of the 
plate Q differs from the above analysis of the single-point case, only in the 
first step. In the case of a multi-point stress we have to construct of ele- 
ments {gQ{x,ar,i), gl{x,ar,i), g2{^:0,r,i)}^=i a basis in the larger deficiency 
subspace Ni, dim Ni = 3m. Due to the presence of singularities of different 
types at different points, the deficiency elements are hnearly independent. 



Concluding remarks on the fitting 
of the model 



The pair of data (21 ) may be used in two different ways: either for calculation 
of the shift of the frequency of SGO and the corresponding perturbation of 
the eigenfunction, under the point-wise stress characterized by the matrix 
M, or, vice versa, for recovering of the data on the localization, the type and 
the intensity of the stress from instrumental observations. 

Indeed, if the geological structure of the tectonic plates at the active zones, 
encoded in matrices Ms attached to the zones, are known, then, theoretically, 
we are able to calculate the eigenf unctions and the eigenfrequencies of the 
plates, taking into account the stress caused by collisions. We are also able, 
theoretically, to construct the deficiency elements for all active zones. Then 
the self-adjoint extension of the bi-harmonic operator on the plate, with the 
point-wise boundary stress, can be constructed, with the corresponding ma- 
trices Mg. The obtained theoretical results can be compared with the results 
of the instrumental measurements. This permits to recover the matrices Mg, 
which characterize the stressed points ai, 02, .... 

Assume that the structure of the plates at the collision point a remains 
unchanged, but the tension is growing linearly with time r as: M(r) = mr. 
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with a matrix coefficient m : N ^ N. Then the formulae ( 20p!9 ), for small 
r, define the derivatives of Xm, (fM with respect to r at the moment r = 0: 



dr T=o 



(Ao - i)K^imoP+(po{xs) 



dr r=o 



Comparing this result with ratios 




A^ 



•mr 



Ao 



r 



r=0 



r 



r=0 



measured experimentally for the amplitude and frequency of SGO, we are 
able to find fit m, and calculate the increment of M 6M = rm. 

Practical experience in analytic perturbations shows, that minor pertur- 
bations affect rather the eigenvalues, than the the shapes of the eigenfunc- 
tions of the spectral problem. Based on this observation we can estimate the 
speed of accumulation of elastic energy S, under the point-wise boundary 
stress depending on the speed of the shift of the eigenvalues (eigenfrequen- 
cies) of SGO and initial distribution of the elastic energy on the modes (/9q 
defined by the corresponding Fourier coefficients (m,v9q): 



with the summation extended only on the eigen-modes which correspond to 
the varying eigenvalues. If only one active zone is involved at a time, then 
only one matrix Mg has to be taken into account, so that the risk of the 
powerful earthquake may be estimated based on the magnitude of 6s = rrris- 
If there are several active zones at the points ai,a2, ■■■ an on the boundary 
of the plate, then the corresponding matrices Mi, M2, . . . M„ can be fitted 
based on observations of SGO in the remote zone during preceding earth- 
quakes which occurred at CLl, CL2, ■ ■ ■ dn- 

Variations of the frequencies and 
the shape of SGO may arise from the stress at any active zone, but usually 
only one active zone is involved at a time. Once the matrices Mi, M2, . . . are 
known, then comparison of the perturbation of the frequencies and the shapes 
of seismo-gravitational modes, SGM (or, probably, seismo-gravitational pul- 
sations, SGP) at the given groups of points in the remote zone with results of 
previous measurements at these points, would allow us to identify the active 
zone where the stress is applied. We presume that the above model gives 
a chance to introduce a useful system into the scope of the experimental 



dS 

dr 



s s 
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data on the seismo-gravitational oscillations in remote zone and use them for 
estimating of risk and localization of the powerful earthquakes. This opens 
an alternative to the statistical methods, see [M], of estimation of risk of 
powerful earthquakes. 

Fitting of the proposed model in reality requires both extended comput- 
ing and a major experimental data base. Because the wavelengths of the 
standing waves, circa 5000 km, on the essential range of periods, 20 min- 1 
hour, dominate the size of the active zones, one may assume, that the straight- 
forward computing with averaged and smoothed data for Young's modulus 
and geometric characteristics of the plates will enables us to obtain a realis- 
tic approximation of the deficiency elements, with singularities at the active 
zones, and to construct the perturbed eigenfunctions of the plates, which 
correspond to SGO. 

More accurate theory requires taking into account realistic boundary con- 
ditions, and the exchange of energy with the liquid underlay and neighboring 
plates. In particular, arising new modes in the spectrum of SGO of the plate, 
transferred, due to the tight contact in active zone, from the neighboring 
plate may be considered as another possible precursor of a powerful earth- 
quake. The choice of realistic boundary conditions has to be done based on 
experimental data interpreted within an appropriate extension of the scheme 
proposed above, with Lame and hydro-dynamical equations involved. We 
postpone discussion of these interesting questions to oncoming publications. 
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